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(54) Abstract Title 

Method for electromagnetic wavefield resolution 

(57) A method of determining the nature of strata beneath the seabed comprises applying an em wavefield to 
the strata and detecting an em wavefield response, and processing the detected response wavefield. The 
processing Is done by resolving the detected wavefield Into upgoing and downgomg components. The 
downgoing component represents reflections from the sea surface which causes interference, while the 
upgoing component represents reflections and refractions from subsea strata. The upgotng component is 
further analysed to derive the nature of the strata. The wavefield frequency is typically between 0.01 - 20 Hz. 
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System and Method for Electromagnetic 
Wavefield Resolution 

The present invention is concerned with electromagnetic data acquisition and 
processing. In particular, the invention relates to a system and method for 
electromagnetic wavefield resolution. 

Marine electromagnetic exploration is an important tool for locating off-shore 
hydrocarbon reserves and monitoring hydrocarbon production for reservoir 
management. One known procedure for marine electromagnetic involves the use of an 
electromagnetic source and receiver cables as described in the present applicants* WO 
01/57555* Electromagnetic energy generated by the source propagates both upwards 
into the water column and downwards through the earth. The downward propagating 
waves are partially reflected and refracted by subsurface layers. The reflected and 
re&acted energy travels upwardly from the subsurface layers and is detected by the 
receiver array. In particular, hydrocarbon filled reservoirs are known to give strongly 
refracted energy which is of interest for hydrocarbon imagmg. 

Electromagnetic exploration however is complicated by waves received at the receiver 
array as downward-travelling reflections and refractions after reflecting and refracting 
off the air/water boundary at the surface. The air/water boundary is an efficient 
reflector and refractor, and thus the waves travelling downwards are difficult to 
differentiate from the upgoing waves from the subsurface. The downward-travelling 
energy is caused both by energy propagating directly firom the electromagnetic source 
to the airAvater boundary and by energy fix>m the subsurface travelling to the air/water 
boundary. 

Reflections and refractions from the sea surface thus are a severe problem. If the sea 
surface reflections and refractions are not properly attenuated, they may interf^e and 
overlap with primary reflections and refractions from the subsurface. 
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It is an object of the present invention to provide a method of processing an EM 
wavefield which minimises this difficulty. 

5 According to the invention a method of processing an electromagnetic (EM) wavefield 
comprises resolving (or decomposing) the wavefield into upgoing and downgoing 
compounds then analysing the upgoing component. Optimal processing, analysis and 
interpretation of electromagnetic data ideally require fiill infonnation about the 
wavefield so that the wavefield can be resolved into its upgoing and downgoing 
10 constituents. 

At a position just above or below the seabed the sea surface reflections and refi-actions 
are always downgoing wavemodes. The reflections and refi:actions of interest fi^om 
the subsurface, however^ are upgoing wavemodes. Resolution (or decomposition) of 
15 the electromagnetic wavefield into upgoing and downgoing constituents just above or 
below the seabed puts the sea surface reflections and reflections into the downgoing 
component whereas the subsurface reflections and refi-actions are contained in the 
upgoing component 

20 Thus, it is a further object of the present invention to provide a technique that resolves 
(or decomposes) the electromagnetic wavefield recorded along one or several receiver 
arrays into upgoing and downgoing wave components. 

Preferably, therefore the wavefield is resolved using the Maxwell Equations: 

25 

AxE(x,t) = ^(z)6tH(x,t) (1) 
A X H (x, t)= [a(z) + e (2) 8,] E (x, t) (2) 

for electric and magnetic fields respectively in an isotropic medium, where; 
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X = (xi, X2, X3) denotes a fixed co-ordinate system with a depth axis positively 
downwards and X3 " 2; n is magnetic permeability, € is magnetic permittivity and a is 
electrical conductivity, whereby \x^Vi (z), e = € (z) and a — u (z); E is the electric field, 
and H is the magnetic field. 

5 

The technique can be used on electromagnetic data recorded on an areal grid or on 
data recorded along a profile (line) or on data recorded as single receiver stations. 
Each recorded component of the electromagnetic wavefieid should be properly 
calibrated before the resolution technique is applied. The calibration ensures that the 
1 0 components of the electromagnetic field satisfy as closely as possible MaxwelFs 
Equations, Preferably, Maxwell's Equations (1) and (2) are transformed using a 
Fourrier transfonn fimction with respect to time and horizontal spatial co-ordinates. 

It is a fiirth^ object of the present invention to provide an approximate technique that 
15 lends itself to adoption in the case of recorded electromagnetic data from individual 
receiver stations, (that is, no summation or integration over receiver stations is 
required). 

Preferably, the upgoing component of the EM wavefieid is derived using the following 
20 formulae: 

U<^^> 1/2 (El - 1/CE H2) (50) 
U<^> = 1/2 (E, - 1/CE H2) (50) 

25 Where U^^^ is the upgoing compound of Ei and Ei is the electric field in a first 

horizontal direction; is the upgoing component of E2 and E2 is the electric field in 
a second horizontal direction; Hi and H2 are the magnetic fields in the first and second 
directions; C is the speed of wave propagation; and E is the complex peraiittivity. 
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Thus, by using Maxwell's Equations, a new method is provided for resolving a marine 



electromagnetic wavefield into upgoing and downgoing wave constituents. The 
effects of the air/water surface can be removed or attenuated through the up/down 
resolution step. The analysis results in expressions where slowness (or wavenumber) 
dependent filters are multiplied with the electromagnetic data Founier transformed to 
a slowness (or wavenumber) domain. After wavefield resolution, the filtered data are 
inverse Founier transformed to a space domain for possible further processing, 
analysis or interpretation* For vertically travelling plane waves the resolution filters 
are independent of slowness: the filters become simple scalers. Wavefield resolution 
then can be carried out directly in a space domain. In this case, the up- and 
downgoing separation is perfomied for each receiver station in the electromagnetic' 
experiment Furthermore, these scalers can be used to resolve approximately the 
electromagnetic wavefield mto upgoing and downgoing components even for the non- 
veitically travelling electromagnetic wavefield. 

For up- and downgoing separation just above the seabed, the resolution filters depend 
on the material parameters of water. For up- and downgoing separation just below the 
sea floor the resolution filters require knowledge of or an estimate of the complex 
wave speed and the complex permittivity (or resistivity, the reciprocal of permittivity) 
of the sea floor material. 

The invention also extends to a method of deteraiining the nature of strata beneath the 
seabed which comprises: applying an electromagnetic (EM) wavefield to the strata; 
detecting an EM wavefield response; and processing the wavefield as described above; 
the nature of the strata being derived Srom the analysis of the upgoing component of 
the detected wavefield response. 

Preferably, the EM field is applied by means of a transmitter located at or near the 
seabed* and the wavefield response is detected by means of a receiver located at or 
near the seabed. Preferably the EM wavefield is transmitted at a fi-equency between 
0.01 and 20Hz. 
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Preferably the transmitter and receiver are dipole antennae through other forms of 
transmitters and receivers can be used. Preferably, the EM wavefield is applied for a 
time in the range 3 seconds to 60 minutes. 

5 

The magnetic measurement nccessaxy may be taken using known magnetotelluric 
instruments. Alternatively, integrated measuring instruments can be used which 
record both magnetic and electric fields. 

10 While the description in this specification mentions the sea and sea bed, it is to be 
understood that these terms are intended to include inland marine systems such as 
. lakes, river deltas etc. 

The invention may be carried into practice in various ways and one approach to the 
15 resolution of the wavefield will now be described in detail^ by way of example, in 
order to illustrate the derivation of formulae the upgoing wavefield component. 

First Maxwell's Equations will be reviewed. Then it will be shown how the 
electromagnetic wavefield can be resolved (or decomposed) into ixpgoing and 
20 downgoing waves. 

A list of the most fi-equently used symbols is given in Appendix A 
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MAXWELL'S EQUATIONS 



We first show how Maacwell's equations can be transfoimed to the frequency- 
horizontal wavenumber domain* Let x — (xi^x^^x^) denote a fixed coordinate system 
with the depth axis positive downwards. For nctationaJ convenience, we will also use 
X3 = z. On the sea floor, assume the material parameters magnetic permeability /i 
and permittivity t as well as the electrical conductivity a do not vary laterally so 
that 

ti^fi(z) ; e=^e(z) ; a = a{z) 

Maxwell's equations for the electric and magentic fields, in conjunction with the 
constitutive relations, for an isotropic medium are given as 

V X E(x,t) = -fi{z)dtHix,t) (1) 
V X H{x,t) = [aiz) £iz)dt] E{x, t) (2) 

where B is electric field, and H is the magnetic field. Introduce the Fourier transform 
with respect to time and horizontal spatial coordinates 

/ / dxdydt txpl-iikixi + k2X2 - £^i)]p(xi ,X2,t) (3) 

with inverse 

2 roo roo roo 

g(xi,X2, t) = -^2^ 7-00 /-ooi ^df^idk2dijjexp[i(k'iXi + ^2X2 - a;t)]G{fci, uj) (4) 
The Fourier transform of equations (1) and (2) gives 




where Ei = Ei(kuk2jZ,cj) is the transformed electric field, etc. In equations (5) to 
(8) we have introduced the complex permittivity 



and 



V we/ 



Pi = ki/u) ; i = l,2 



(9) 



(10) 



Matrix vector diSerential equation 



. Equations (5) to (8) can be written as an ordinary maitrix-vector differential equa- 



tion 

where the wave vector b is a 4x1 column vector 

El 



(11) 



6 = 



(12) 



E2 
Hi 

aad the system matrix .4 is a 4x4 matrix partitioned into four 2x2 submatrices of 
which the diagonal ones are zero, 

0 Ai 
A2 0 

The submatrices Ai and A2 are synunetric 



A = 



(13) 



Ai = 



a -Si. _Ei£i 



; A2 = 



M J 



(14) 



At and A2 are functions of the parameters in Maxwell's equations (and therefore, 
functions of z) and of p.- 
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DECOMPOSITION INTO UP- AND DOWNGOING WAVES 

For the decomposition of the electromagnetic field into up- and downgoing waves 
it is necessary to find the eigenvalues and eigenvectors of the system matrix A for 
^ven wavenumbers and frequencies. The wave vector b can be decomposed into up- 
and downgoing waves 

t/; = [l7^Z?T. (15) 

where 17^ = [Uu U2] and = [Ui, i?2l, by the linear transformation 

b ^ Lw, (16) 

where L is the locaJ eigenvector matrix of A. (i.e., each column of L is an eigenvector). 
Since L is the eigevector matrix of A it follows that 

where A is the diagonal matrix of the corresponding eigenvalues of A: 

A = diag[— Ai , — A2 , Ai , A2] (1 7) 

ESg^values of A 

The eigenvalues of A are 

^,=A2 = 9=(c-2-p=)^/=' (18) 

where 

c-= = c/i (19) 



(20) 
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Eigenvector matrix of A 

The eigenvector matric of A can be given as 



L = 



■with inverse 













eg 








_ Vim 










<« 




0 


-1 


0 


-1 


1 


0 


1 


0 







0 


1 






-1 


0 






0 


1 






-1 


0 



(21) 



(22) 



Upgoing and downgoing waves 

From equation (16) upgoing and downgoing -waves are given by 



that is 



^> = 2 [ — 



21 q - q 
As is shown below, I7i, £^2, and Z>2 have been defined such that 



(23) 



(24) 
(25) 
(26) 
(27) 



(28) 



This implies that Ui and 2?i are the upgoing and dofsragping constituents of Hu 
respectively, whereas U2 and D2 are the upgoing and downgoing constituents of H2, 



10 

respectively- The scaling of upgoing and downgoing waves is however not unique. 
We will show below that the upgoing and downgoing waves defined in equation (27) 
can be scaled such their sum yields upgoing and downgoing constituents of the fields 
El and £2- The upgoing constituents of Hj^ H^^ E\ and £?2 will not contain the 
downgoing reflections and refractions caused by the sea surface- After decomposing 
the measured electromagnetic field into upgoing and downgoing wave fields the sea 
surface reflections and refi-actions will belong to the downgoing part of the fields- The 
upgoing and downgoing wavefields are inverse Fourier transformed to space domain 
using equation (4). 

UPGOING AND DOWNGOING CONSTITUENTS OF Hi AND Hq^ 

Equation (28) is easily verified hy summation of Ui and A , and V2 and as given 
in equation (27). Therefore, the wave fields U\ and Di are interpreted as upgoing and 
downgoing constituents of the magnetic field component Hu whereas the wave fields 
U2 and D2 are interpreted as upgoing and downgoing constituents of the magnetic 
field component H2* We introduce the notation 

J7<^*> = £/i ; = A = ffi - 17<^^) (29) 

U^^'^^^U^ ; D^^^^^D2 = H2-U^''^^ (30) 

so that 

J?i = t/t^i)4-X><'^»> ; J?2 = i7f^O + I?f^=^5 (31) 
In particular, the upgoing constituents (see equations (24) and (25)) are of interest 
i7(^>) = i [^1 - y (piP2^i - ^£2)] (32) 

C^^^^ = ^ [^2 -f- y (P1P2£2 - 9^£?:)] (33) 

Equations (32) and (33) are the most general formulas for electromagnetic wavefield 
decomposition of the magnetic field components into upgoing waves* The schemes 
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require the receiver stations to be distributed over an areal on the sea bed so that the 
dectromagnetic Tvavefield can be tranformed to the slowness domain. The decompo- 
sition schemes (32) and (33) are valid for a 3D inhomogeneous earth. 

Special case: p2 = 0 

When the electromagnetic experiment is run along a single profile electromagnetic 
data are available along a line only. The magnetic field components and ff^ then 
can be properly decomposed into its upgoing and downgoing waves under the 2,5D 
earth assupaption (no variations in the medium parameters of the earth in the cross- 
profile direction). Without loss of generality, orient the coordinate system so that 
the electromagnetic wavefield propagates in the xj^xs-plane such that p2 = 0* Then, 
52 = c"^ g = ft» inserted into equation (32) gives 

C^(HO = l(fr, + c^j,e£i) (34) 

Equation (34) shows that to remove the downgoing reflected and refracted energy 
from the Hi magnetic field it is necessary to combine the Hi recording with a scaled 
(filtered) E2 electric field recording. Similarly, the upgoing component of the H2 field 
is 

Equations (34) and (35) are strictly valid under the 2.5D earth assumption. However, 
for single profile data over a 3D earth equations (34) and (35) still can be used as 
approximate methods to attenuate the downgoing energy on the magnetic H^ and H2 
components. 

Special case: pi = P2 =^ 0 

The special case of vertically traveling electromagnetic plane waves with pi = 
P2 = 0 such that = q2 ^ g — c"* yields by substitution into equations (32) and 
(33) 
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i7(^v)wi(jyi+c€^) (36) 
J7<^2) = i(H2_C££;x) (37) 

Even though equations (36) and (37) are strictly valid only for vertically traveling 
plane waves as a decomposition method for the magnetic components, they can be 
a useful approximation for Tvavefield decomposition also for non-vertically traveling 
plane waves as well as for the full magnetic Hi and H^, fields. Note that since the 
scaling factor applied to the electric components does not depend on slowness, equa* 
tions (36) and (37) can be implemented in space domain. In this special case, Hi or 
H2 magnetic data recorded on each receiver station are processed independently 

UPGOENG AND DOWNGOING CONSTITUENTS OF AND Ik 

By properly scaling the upgoing and downgoing waves Ui^ C/2> Di and D2 we can 
find the upgoing and downgoing constituents of the fields E\ and £2. The scaling 
must chosen to give 

El = i7^^>) + i)(^^5 (38) 

(39) 



with 



C^{*^>) = t;p> + t/pi) (40) 

£><^^> = Di^>> + i?f^> (41) 

U{B2) = lf[E*) _^ yiB^) (42) 

Z)(^) = Z)(^>+2?<^> (43) 



Introducing 



CK^O^^C/, i7P^J = -^t/2 ; Df'^^-^Dr ; P^^^' = (44) 
€? tq eg eq 

we find that equation (38) is fulfilled^ and that 

Ui^) = C^P> + C/i^'J = l[^i + ^(PiP2^i - (45) 




eq eq eg 




D2 (46) 



we find that equation (39) is fulfilled, and that 



(47) 



Equations (45) and (47) are the most general formulas for electromagnetic wavefield 
decomposition of the electric field components into upgoing ^ves. The schemes re^ 
quire the receiver stations to be distributed over an areal on the sea bed so that the 
electromagnetic wavefield can be tranformed to the slowness domain. The decompo- 
sition schemes (45) and (47) are valid for a 3D inhomogeneous earth. 

Special case: P2^0 

When the electroznagnetic experiment is run along a single profile electromagnetic 
data are available along a line only. The electric field components and E2 then 
can be properly decomposed into its upgoing and downgoing waves under the 2.5D 
earth assumption (no variations in the medium parameters of the earth in the cross- 
profile direction)* Without loss of generality, orient the coordinate system so that 
the electromagnetic wavefield propagates in the a;i,3:3-plane such that P2 = 0. Then, 
92 = c"*, g ft , inserted into equation (45) gives 



Equation (48) shows that to remove the downgoing reflected and refracted energy 
from the £^1 electric field it is necessary to combine the Ei recording with a scaled 



Equations (48) and (49) are strictly valid under the 2.5D earth assumption. However, 
for single profile data over a 3D earth equations (48) and (49) still can be used as an 
approximate method to attenuate the downgoing energy on the electric £1 and B2 
components. 




(48) 



(filtered) H2 magnetic field. Similarly, the upgoing component of the B2 field is 




(49) 



Special case: j?i =: p2 = 0 
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The special case of vertically tiaveling electromagnetic plane waves with pi = 
Pa = 0 such that 9i = ^2 = 9 = c"^ yields by substitution into equations (45) and 



Even though equations (50) and (51) are strictly valid only for vertically traveling 
plane waves as a decomposition method for the electric components, it can be a 
useful approximation for wavefield decomposition also for non-vertically traveling 
plane waves as well as for the fall electric JSi and fields. Note that since the 
scaling factor applied to the magnetic components does not depend on slowness, 
equation (50) can be implemented in space domain. In this special case, or £2 
electric data recorded on each receiver station are processed independently. 



(47) 




(50) 




(51) 
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Appendix A 




A: 


system matrix 


b: 


^ve vector cantaining electromagnetic fields 


tot . 


wave vector containing upgoiDg and downgoing 


L: 


eigenvector matrix of A 


B: 


magnetic flux denaty 


Hz 


magnetic field; H {Huffttffz) 


D: 


eiectiic displacement field 


E: 


electric field; B = {^i,£i,JE^) 


Jz 


current deosily 



SB — (xi.xztXa): Cartesian coordinate 

mm 



cz 



k: 

kjZ 

fe: 
Pi: 

P2: 

P- 

9- 

z: 



upgoing component of £1; £!| == + jj{Bt) 
downgoing compaoent of Ea 
upgoing component of ^ = Z7(«>) + x)(£i) 
downgoing component of B2 
upgoing component of iTi; fli = £/^(fft) + 
downgoing component of Hi 
upgoing component of ffi = U^^*) + D{ff») 
downgoing component of £2 
Speed of wave propagation; c = (/je)-^ 
Wavennmber; k = u/c 
Horizontal wavenumber conjugate to Xi 
Horizontal wavenumber coiuugate to 
BoxizontaJ slowness pi s ki/ut 
Horizontal slowness p2 «= kz/u 

Vertical slowness; g = yj<r^—j^-i^ 
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p^i volume electric charge density 

p: resistivity; the reciprocal of resistivity is conductivity 
ei permittivity 

c: complex permittivity, c = e (l + ^) 
lit magnetic permeability 

£T: electrical conductivity; the reciprocal of conductivity is resistivity 

A: eigenvector 

u: circular frequency 

dti temporal derivative; 9t = ^ 

d}i spatial derivative; dx = 

atj: spatial derivative; ^ - ^ 

ds: spatial derivative; ^ = ^ 
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Claims 

L ■ A method ofdetenninmg the nature ofstrata beneath the seabed wHch 
comprises: applying an electromagnetic (EM) wavefield to the strata; detecting an EM 
wavefield response; and processing the wavefield by resolving the wavefield to 
produce upgoing and downgoing components then analysing the upgoing component, 
thereby deriving the nature of the strata. 

2. A method as claimed in Claim 1, in which the EM field is applied by means of 
a transmitter located at or near the seabed. 

3. A method as claimed in Claim 1 or Claim 2, in which the EM wavefield is 
transmitted at a ficquency between 0.01 and 20Hz. 

. 4. A method as claimed in any preceding Claim, in which the EM wavefield is 
transmitted at a wavelength between IS and SOS. where S is the thickness of the 
overburden above the considered strata. 

5. A method as claimed in any preceding Claim, in which the transmitter is a 
dipole antenna. 

6. A method as claimed in any preceding Claim, in which the wavefield response 
is detected by means of a receiver located at or near the seabed. 

7. A method as claimed in any preceding Claim, in which the receiver is a dipole 
antenna. 

8. A method as claimed in any of Claims 1 to 6, in which the receiver comprises a 
detector pair consisting of means for detecting the electromagnetic wavefield and 
means for detecting a magnetic field- 
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9, A method as claimed in Claim 8 in which the detector pair is housed in a single 
unit* 

10. A method as claimed in Claim 8 or Claim 9, in which data from one of the 
detector pairs is used to resolve the field detected by the other of the detector pairs. 

1 L A method as claimed in any preceding Claim, in which the wavefield is 
detected using a plurality of receivers airanged over an area of the seabed. 

12. A method as claimed in Claim 11, in which the receivers are arranged in a line. 

13. A method as claimed in Claim 11 or Claim 12, in which data from an array of 
receivers is used to resolve the wavefield. 

14. A method as claimed in Claim 1 1 or Claim 12, in which data from each 
receiver independently is used to resolve the wavefield. 

15. A method as claimed in any of Claims 6 to 14, in which the or each receiver is 
moved to a different location while the transmitter is maintained stationary, and a 
further EM wavefield is applied, detected and processed. 

16. A method as claimed in any preceding Claim, in which the wavefield is 
resolved using the 2.SD assumption that there are no variations in the medium 
parameters in the cross-profile direction. 

17. A method as claimed in any preceding Claim, in which the resolution and 
analysis are carried out on the assimiption of vertically travelling plane waves. 



